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ABSTRACT 

We present the results of ./V-body simulations in which we take the masses, posi- 
tions and velociti es of sink particles from five pairs of hydrodynamical simulations of 
star formation bv lDale et al.l ()2012l . [20131 ) and evolve them for a further lOMyr. We 
compare the dynamical evolution of star clusters that formed under the influence of 
mass-loss driven by photoionization feedback, to the evolution of clusters that formed 
without feedback. We remove any remaining gas and follow the evolution of structure 
in the clusters (measured by the Q-parameter) , half- mass radius, central density, sur- 
face density and the fraction of bound stars. There is little discernible difference in 
the evolution of clusters that formed with feedback compared to those that formed 
without. The only clear trend is that all clusters which form without feedback in the 
hydrodynamical simulations lose any initial structure over lOMyr, whereas some of 
the clusters which form with feedback retain structure for the duration of the sub- 
sequent iV-body simulation. This is due to low e r init i al den sities (and hence longer 
relaxation times) in the clusters from lDale et al.l (|2012l [2013h which formed with feed- 
back, which prevents dynamical mixing from erasing substructure. However, several 
other conditions (such as supervirial initial velocities) also preserve substructure, so 
at a given epoch one would require knowledge of the initial density and virial state of 
the cluster in order to determine whether star formation in a cluster has been strongly 
influenced by feedback. 
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1 INTRODUCTION 

The conversion of molecular gas to stars results, at 
least initially, in stars g eometrically grouped into clus- 
ters and associations (e.| 



ters and associations (e.g. lLada fc Lada 120031 ; lLadal l2010l : 
iGieles fc Portegies Zwartl |201ll ). The transition from em- 
bedded clusters whose dynamical behaviour is dominated 
by gas, to systems devoid of gas, is still not well under- 
stood. Stellar feedback, particularly in the form of pho- 
toionizing radiation, winds and supernovae from OB-type 
stars, but also photodissociating radiation and jets from low- 
and intermediate-mass stars, are strong candidates for the 
agency of gas expulsion an d much theoretica l work has been 
devoted to their study (e.g. IWhitworth|[l979l: lMatznerll2002i : 
iLi fc Nakamurall2006l : lKrumholz. Klein fc McKeell2012i '). 

Although the dynamics of gas-free clusters are, by defi- 
nition, governed by gravity (including tidal interactions with 
other bodies) and stellar evolution, the prior hydrodynamic 
phase and the transition from that phase to a purely grav- 
itational system are of crucial importance, since together 
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they determine the initial stellar masses, positions and ve- 
locities for the subsequent dynamics-dominated evolution of 
the cluster. 

In a re cent series of papers, Dale. Ercolano fc Bonnelll 

(|2012h and iDale. Ercolano fc Bonnelll (|2013l ) modelled the 
formation of massive clusters with and without the influence 
of feedback. These authors studied the influence of mass-loss 
driven by photoionizing radiation from massive stars (m > 
20 Mq ) on the dynamical evolution of their model clouds by 
comparing each simulation with a corresponding control run 
in which feedback from the massive stars was switched off. 
The simulations were evolved for ~3 Myr, during which time 
the clusters generally remained substructured, although dif- 
ferences between the runs with and without feedback were 
generally limited to the respective initial mass functions, 
which were somewhat more top-heavy for the simulations 
where star formation was not regulated by feedback. 

Recently, iV-body simulations of cluster evolution have 
been used to demonstrate that the initial conditions of star 
formation can in some cases be discerned from a process 
of 'reverse engineering' - comparing snapshots from the 
simulations to observed clusters and re-calibrating the ini- 
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tial conditions of the simulations. By assuming realistic ini- 
tial morphologies for clusters (i .e. clumpy, or substruct ured : 
ICartwright fc Whitworth1l2004ISanchez fc Alfardl2009h an d 



2 INITIAL CONDITIONS 

2.1 Embedded cluster formation 



subvirial velocities ( Peretto et alj 20061 ; Furesz et al.l 2008), 
it is possible to dynamically mass-segregate clusters on very 
short timescales (< 1 Myr) - the more clumpy the initial 
conditions, the mo re likely mass segregation is to occur 
jAllison et alj|20ld). 

Additionally, lAllisonl (|2012T) showed that the more sub- 
virial and substructured a cluster is, the more stars are 
ejected from the cluster, and the larger their velocities are. 
The initial virial state of a cluster c an also be inferred from 
the amount of substructure present. iParker fc Meverl (|2012T ) 
demonstrate that a subvirial, substructured cluster loses its 
substructure on very short timescales, whereas a supervirial 
(expanding) cluster retains substructure in most cases. 

In principle, it may be possible to find differences be- 
tween clusters that form with feedback and those that form 
without by examining their long-term dynamical evolution. 
In this paper, we perform A-body simulations in which we 
take as initial conditions the end results of SP H hydrody- 
namical simulations of embedded clusters from Dale et al.l 
(|2012h and lDale et all (|2013l ). Here, we have chosen the five 
pairs of feedback and control simulations from these works 
in which individual stars could be resolved. We assert that 
the remaining gas is expelled instantaneously by supernovae, 
and treat the resulting clusters as A-body problems. 

Earlier work on the impact of instantaneous gas ex- 
pulsion on clusters suggested that on the removal of back- 
ground gas po tential, clusters quickly become unbound and 



ground gas potential, clusters quickly become unbound and 
disperse (e.g. [Tutukov 19781; lLada. Margulis fc Dearborn! 



11984 lGoodwmlll997l :l Bastian fc Goodwin! I2006T 1 . However , 
if a cluster is initially substructured, then the subclus- 
ters (which contain the majority of stars) can be gas-poor 
and the removal of the potential has onl y a minimal ef- 
fect o n the global dynamics of the cluster l|Kruiissen et al.l 
12012]) . This is due to the stars a nd g as being decoupled 
l Offner, Hansen fc Krumholzll2009D . and lSmith etafl (|201lh 
find that the distribution of initial stellar positions and ve- 
locities determines whether the star cluster survives, rather 
than the star formation efficiencjQ. 

Here, we remove the gas before the A-body evolution 
and rather than examine the effect of removing the gas po- 
tential, we determine what influence (if any) the previous 
gradual gas expulsion by photoionization has on the sub- 
sequent (gas free) behaviour of the systems. We briefly de- 
scribe our initial conditions and numerical methods in Sec- 
tion 2. The results of the A-body calculations follow in Sec- 
tion 3 and our discussion and conclusions follow in Sections 
4 and 5. 



1 Another mechanism to limit the effects of gas expulsion is a 
density— driven star formation efficiency g radient in the cluster 
iAdamsll200(i iParmentier fc Pfalznenl2013h , although this mech- 
anism assumes spherical symmetry, which would require clusters 
to have lost any initial substructure. 



In a re cent series of papers. fPale et all l|2012T ) and lDale et al.l 
(|2013h presented a series of Smoothed Particle Hydrodynam- 
ics (SPH) simulations of the evolution of embedded clusters 
in GMCs with two initial cloud virial ratios of 0.7 and 2.3 
(i.e. bound or unbound), and a range of masses (10 4 M o - 
10 6 Mq) and radii between 2.5 and 100 pc. The initial den- 
sity profiles of the gas clouds are all close to p oc r~ 2 , since 
the equation of state used is close to being isothermal. We 
concentrate here on the lower-mass (l- 3xl0 4 M ) clouds 
from iDale et all (|2012T ) and iDale et alj (12013T ). since these 
are the simulations in which individual stars could be re- 
solved. 

IDale et alj (|2012l . l2013h allowed their model clouds to 
form a few ionizing sources, then split each simulation in 
two. In one version, feedback from these massive stars is 
switched on, while in the second (control) run, feedback is 
left switched off. The formation and evolution of the clusters 
is then followed for ~3 Myr, in an attempt to understand the 
isolated effect of this form of feedback before the detonation 
of any supernovae. 

Before the onset of feedback, IDale et alj (|2012l . 120131 ') 
examined the clusters' st r uctur e using the Q-para meter 
jCartwright fc Whitworthl |2004| ; ICartwrightl 120091 ) and 
found that all the clusters have some initial substructure, 
which is retained at the end of the SPH simulations (some 
clusters become more substructured, others less substruc- 
tured). TJri£i££c^iMstent ndth the simulations of star forma- 
tion bv lGirichidis et alj |201ll ), in which the majority of the 
clusters formed in their simulations (from a range of initial 
cloud density profiles) were substructured dGirichidi s et alj 
|2012| . see their table 3). 

In the simulations bv lDale et all (|2012l . |20lj ), the dif- 
ferences between the runs with and without feedback were 
generally limited to the respective initial mass functions, 
which were somewhat more top-heavy for the feedback-free 
simulations. The clouds lost a wide range of fractions of 
their mass during the formation of stars - from negligible 
amounts to more than half - and mass-loss had a largely 
commensurate influence on the structure and dynamics of 
the embedded clusters that formed. 

In the work we present h ere, we only use the simulations 
from [Dale et al.l(|2012l . l2013l ) in which individual stars could 
be resolved (down to the lowest possible stellar mass allowed 
by the resolution limit of the simulation - 0.5 Mq). This 
gives us a total of five sets of simulations which we will use as 
the initial conditions of subsequent A-body simulations, to 
compare the effects of dynamical evolution between clusters 
that formed under the influence of feedback, and those that 
formed without. In Table[T]we summarise these five pairs of 
SPH simula tions, and present the Run ID from IDale et alj 
l|2012l . I2013T I - the prefix 'IP indicates an unbound cloud - 
as well as the initial physical parameters of the cloud, the 
number and total mass of stars formed, and the initial and 
final Q-parameters in the SPH simulations. 



2.2 A-body evolution 

We now take the final states of five pairs of simulations 



We now take tne hna i stat es or nve pairs 1 
from IDale et alj (|2012r i and IDale et all (j2013T ) 



and assume 
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Table 1. A summary of the five different pairs of smoothed particle hydrodynamics (SPH) simulations used as th e input initial cond itions 
of our TV-body integ rations. The values in the columns are: the simulation number, corresponding Run ID from [Dale et al, | i2012l , D12) 
or iDale et al] l|2013l . D13), status of feedback in the SPH simulation ('on' or 'off'), the paper reference, the initial virial ratio of the 
original clouds oj^J 1 (to distinguish bound from unbound clouds), the initial radius of the cloud in the SPH simulation (il c i ou d), the 
initial mass of the cloud (-Mcloud)i the number of stars that have formed at the end of the SPH simulation (Altars), the mass of this 
cluster (M c i uster ), the Q-parameter in the SPH simulation at the time feedback is initiated in the feedback runs (Sf™), and the final 
Q-parameter in the SPH simulations (<2|^ H ). 



Sim. No. 


Run ID 


Feedback 


Ref. 


"init 


^•cloud 


A^cloud 


-/Vst ars 


-^cluster 


qSph 

^init 


QSPH 


1(a) 


J 


On 


D12 


0.7 


5 pc 


10 000 M Q 


685 


2205 M Q 


0.53 


0.60 


1(b) 


.1 


Off 


D12 


0.7 


5pc 


10 000 M Q 


578 


3207 Mq 


0.53 


0.49 


2(a) 


I 


On 


D12 


0.7 


10 pc 


10 000 M 


168 


805 Mq 


0.42 


0.38 


2(b) 


I 


Off 


D12 


0.7 


10 pc 


10 000 Mq 


186 


1270 M 


0.42 


0.72 


3(a) 


UF 


On 


D13 


2.3 


10 pc 


30 000 M 


76 


836 M q 


0.59 


0.55 


3(b) 


UF 


Off 


D13 


2.3 


10 pc 


30 000 M 


66 


1392 Mq 


0.59 


0.77 


4(a) 


UP 


On 


D13 


2.3 


2.5 pc 


10 000 M Q 


346 


1957 M q 


0.47 


0.57 


4(b) 


UP 


Off 


D13 


2.3 


2.5pc 


10 000 M Q 


340 


2718 Mq 


0.47 


0.49 


5(a) 


UQ 


On 


D13 


2.3 


5 pc 


10 000 Mq 


80 


648 M q 


0.42 


0.46 


5(b) 


UQ 


Off 


D13 


2.3 


5 pc 


10 000 Mq 


48 


723 Mq 


0.42 


0.70 



that the detonation of the first supernova instantaneously 
removes any remaining gas from both the feedback and non- 
feedback simulations, and evolve the resulting gas-free sys- 
tems with an A^-body code. This allows us to gauge the in- 
fluence of gradual and prolonged ionization-driven mass-loss 
on the subsequent long-term evolution of the clusters. 

As discussed in IDale et al l (|2012T ) and IDale et all 
l|2013l ). and Section [5TTJ the lowest mass star which forms in 
the SPH simulations has a mass of ~0.5 Mq. This resolution- 
limited IMF should not affect the results presented here, as 
we are looking for differences between the evolution of sim- 
ulations with the same mass resolution. 

We evolve the gas-less clusters using the 4 th order 
Hermite-sc heme integrator kira within the S tarlab environ- 
ment (e.g. iPortegies Zwart et al.lll999l , 120011 ). We take the 
masses, positions and velocities of the sink-particles from 
the SPH simulations and place these directly into the N- 
body integrator. In the majority of the SPH runs the stars 
are in virial equilibrium, or slightly sub-virial, at the end 
of the simulation (i.e. the initial conditions for the A^-body 
integration). However, the initial conditions for simulations 
UF, UP and UQ were globally unbound, so that one might 
expect the stars and clusters formed to be in an unbound 
configuration. However, some parts of the globally unbound 
clouds become bound due to high velocity gas flows colliding 
and radiating away kinetic energy, which tends to occur in 
the dense areas of the clouds where most of the stars form. 
Therefore, in practice, these unbound clouds form stars that 
are roughly in virial equilibrium (with virial ratios ranging 
from 0.4 - 0.7), apart from Run UF (with feedback), which 
is highly supervirial, with a virial ratio of 1.9. 

In their recent paper, iMoeckel et all |2012T ) used 
the output o f the SPH simulation of star formation by 
iBonnell et ai1l|2003l ,l 2008) as the input initial conditions for 
A^-body simulations to study the subsequent gravitational 
evolution of the sub-clusters which formed in the SPH sim- 
ulation. They deliberately scaled the input velocities of their 
sink-particles to be in virial equilibrium. In this work, we will 



use as a default the original (non-scaled) velocities from the 
SPH simulations, but we also run the simulations with the 
stellar velocities scaled to virial equilibrium in order to make 
a comparison. 

The simulation pairs are then evolved for 10 Myr, with- 
out a background gas potential. 

The simulations contain several stars with masses 
> 20 Mq which are likely to evolve over the 10 Myr 
duration of the A^-body integration. For this reason, we 
use the Se Ba stellar evolution package in the Sta rlab en- 
vironment (|Portegies Zwart fc Verbunq|l996l . l2012h . which 
provides look-up tables for the evolution of stars ac- 
cording to the time dependent mass- radius rela ti ons i n 
lEggleton. Fitchett fc Tout] j 19891 ) and iTout et all l|l99rj) . 
Typically, SeBa updates the evolutionary status of stars on 
shorter timescales than the timestep in the kira integra- 
tor, although for extremely close systems a lag of up to one 
timestep can occur. 

Due to the relatively young end-time of the simulations 
we do not impose an external galactic tidal field on the sim- 
ulations. 



3 RESULTS 

In this Section, we describe the fu rther evolu t ion o f the 
clusters formed in the simulations in IDale et al l l|2012l ) and 
IDale et all l|2013h . using A^-body simulations. We compare 
the evolution of cluster structure (as measured by the Q- 
parameter), mass loss due to stellar evolution, half-mass 
radii, central volume densities, surface densities, the frac- 
tion of bound/unbound stars and search for mass segrega- 
tion. Our results are necessarily only valid f or the initial 
conditions adopted bv lDale et all |2012| , |2013| ). A summary 
of our results is provided in Table [2] 
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Table 2. A s ummary of t he res ult s of our iV-body in tegrations. The values in the columns are: the simulation number and corresponding 
Run ID from lDale et al, I d2012h or lDale et al, | i2013h ('a' corresponds to the SPH run with feedback switched on, and 'b' corresponds to 
the SPH run where feedback was switched off), the initial mass of this cluster before W-body integration (^cluster i)> the final mass after 
lOMyr of ./V-body integration (Mcluster,f ), the initial half-mass radius ()"i/2,i), the final half-mass radius after lOMyr (ri/2,f), initial and 
final volume densities (pi and pf), initial and final median surface densities (Ej and Sf) and the initial and final Q-parameters (Qi and 

Sf)- 



Sim. No. 


-^cluster, i 


^cluster, f 


r l/2,i 


r l/2,f 


Pi 




Pi 




Si 






Si 


S f 


J, 1(a) 
J, 1(b) 


2205 M Q 
3207 M q 


1857 M 
2531 M 


1.69 pc 
1.20 pc 


5.30 pc 
15.8pc 


54 M pc" 
220 M pc" 


3 
-3 


1.5M pc~ 
0.08 M pc" 


3 
-3 


141 stars pc" 2 
4518 stars pc -2 


2 stars pc - ' 
0.4 stars pc~ 


J 

■2 


0.60 
0.49 


1.89 
1.91 


I, 2(a) 
I, 2(b) 


805 Mq 
1271 M 


640 M 
751 M 


3.59 pc 
0.70 pc 


17.6 pc 
8.72 pc 


2.08 M pc 
447 M pc" 


-3 
-3 


0.01 M pc" 
0.14M© pc" 


-3 
-3 


83 stars pc — 2 
102 stars pc — 2 


0.1 stars pc~ 
0.3 stars pc~ 


-2 
■2 


0.38 
0.72 


0.79 
1.39 


UF, 3(a) 
UF, 3(b) 


836 M Q 
1392 M 


511 M 
410 M 


6.85 pc 
1.49 pc 


32.4 pc 

20.5 pc 


0.31 M pc 
50 M pc~ 


-3 
3 


0.002 M pc 
0.006 M pc 


-3 
-3 


0.6 stars pc — 2 
6 stars pc — 2 


0.01 stars pc 
0.01 stars pc 


-2 
-2 


0.55 
0.77 


0.74 
1.01 


UP, 4(a) 
UP, 4(b) 


1957 Mq 
2718 M 


1587 M 
1765 M 


3.85 pc 
3.87 pc 


12.4 pc 
19.6 pc 


4.09 M pc 
5.53 M pc 


-3 
-3 


0.09 M pc" 
0.03 M pc" 


-3 
-3 


24 stars pc — 2 
250 stars pc -2 


0.2 stars pc~ 
0.2 stars pc~ 


■2 
-2 


0.57 
0.49 


1.27 
1.40 


UQ, 5(a) 
UQ, 5(b) 


648 M 
723 M 


485 M 
337 M 


5.13 pc 
1.32 pc 


25.7pc 
13.7pc 


0.57 M pc 
37.3 M pc 


-3 
-3 


0.003 M pc 
0.02 M pc" 


-3 
-3 


2 stars pc — 2 
6 stars pc — 2 


0.04 stars pc 
0.02 stars pc 


-2 
-2 


0.46 
0.70 


0.72 
0.93 



3.1 Cluster morphologies and initial structure 

In Figs.[T]and[2]we show the morphologies of Run I with and 
without feedback, respectively. From top to bottom, we show 
the initial morphology before iV-body evolution (0 Myr) and 
after 1 Myr and 5 Myr. Each panel shows a render of the clus- 
ter, along with projections of the cluster onto the x, y and 
z planes. The simulation with feedback switched on (Fig. [TJ 
exhibits a high degree of substructure, which is still evi- 
dent after 1 Myr (middle panel) and to a lesser degree after 
5 Myr. On the other hand, the simulation without feedback 
(Fig. [2} appears to be slightly substructured initially, before 
assuming a smooth, centrally concentrated morphology. 

Whilst Run I is typical of the differences between all 
pairs of simulations, it is useful to compare the initial mor- 
phologies for the feedback versus no-feedback initial con- 
ditions for all simulations. In order to quantify structure , 
we adopt the Q- parameter (|Cartwright fc Whitworthll2004l ; 
ICartwrightJl2009l ), which compares the mean length of a min- 
imum spanning tree (MST) joining all the stars in a cluster, 
m, to the average separation between stars in the cluster, s: 



In two dimensions, a cluster is substructured if the Q- 
parameter is less than 0.8 (lower values correspond to a 
higher degree of structure), or centrally concentrated if the 
Q-parameter is greater than 0.8 (higher values correspond 
to a more centrally concentrated morphology). In all cases, 
we project our clusters along the z-axis, but the results are 
robust against other projections. 

For Run J, the Q-parameter is 0.60 at the start of the 
iV-body integration for the simulation with feedback, and 
is 0.49 for the simulation with no feedback. On the other 
hand (and as is visible by comparing the top panels of 
Figs. [T] and [2)| , for Run I the simulation with feedback is 
initially more structured, with a Q-parameter of 0.38, as op- 
posed to a Q-parameter of 0.72 for the run with no feedback. 
In the remaining simulations, Run UF has Q-parameters of 
0.55 (with feedback) and 0.77 (no feedback); Run UP has Q- 



parameters of 0.57 (with feedback) and 0.49 (no feedback); 
and Run UQ has Q-parameters of 0.46 (with feedback) and 
0.70 (no feedback). 

In summary, there is no clear trend for clusters with 
feedback to be more or less structured than those without 
at the end of the SPH simulations (and before iV-body in- 
tegration) . 



3.2 Evolution of cluster structure 

We now examine the evolution of the structure of the model 
clusters over the duration of the iV-body integra tion, using 
the Q-parameter l|Cartwright fc Whitworthll2004l ). In Fig. [3] 
we compare the evolution of the Q-parameter with time in 
the five pairs of simulations. The simulations that formed 
with feedback are shown by the solid lines, and the simula- 
tions that formed without feedback are shown by the dashed 
lines. The colours correspond to the following simulations; 
red-Run J, green-Run I, dark blue-Run UF, cyan-Run UP 
and magenta-Run UQ. As detailed in Section 13.11 all clus - 
ters which form in the simulations of lDale et all (|2012l , l2013r ) 
have substructure initially (i.e. Q < 0.8), but there is no 
trend towards a particular regime of Q as a consequence of 
feedback being switched on or off. 

In all cases, the clusters lose substructure as they evolve, 
so that their Q-parameters increase. The clusters result- 
ing from simulations without feedback (the dashed lines in 
Fig. [3]l lose substructure more quickly, and all of these no- 
feedback runs attain a centrally concentrated morphology 
with final Q-parameters of 1.91 (Run J), 1.39 (Run I), 1.01 
(Run UF), 1.40 (Run UP) and 0.93 (Run UQ). 

Conversely, the clusters resulting from simulations with 
feedback lose their structure more slowly, and two simula- 
tions (Run UF and Run UQ) remain substructured even 
after 10 Myr of A^-body integration (the solid blue and ma- 
genta lines in Fig. with final Q-parameters of 0.74 and 
0.72, respectively. Run I with feedback (the solid green line) 
retains substructure for 8.3 Myr, with a final Q-parameter of 
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Run I, feedback ON, time: Myr Run I, feedback OFF, time: Myr 




Run I, feedback ON, time: 1 Myr Run I, feedback OFF, time: 1 Myr 




20 -20 20 -20 



Run I, feedback ON, time: 5 Myr Run I, feedback OFF, time: 5 Myr 




60 -60 



Figure 1. Evolution of cluster morphology for Run I with feed- Figure 2. Evolution of cluster morphology for Run I with no 

back after 0, 1 and 5 Myr of JV-body evolution, with stars colour- feedback after 0, 1 and 5 Myr of TV-body evolution, with stars 

coded by mass. Each panel shows a render of the cluster, along colour-coded by mass. Each panel shows a render of the cluster, 

with projections of the cluster onto the x, y and z planes. The along with projections of the cluster onto the x, y and z planes, 

cluster is still substructured after 5 Myr. The small amount of initial substructure is erased by dynamics 

and the cluster reaches a smooth, centrally concentrated mor- 

© 2013 RAS, MNRAS 000, [HM] phology after 5 Myr. 
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0.8 (i.e. straddling the boundary between substructured and 
centrally concentrated). Runs J and UP with feedback (the 
solid red and cyan lines) lose their substructure within the 
first 2Myr of TV-body integration, with final Q-parameters 
of 1.89 and 1.27, respectively. 

Not only do the clusters resulting from simulations with 
feedback lose structure more slowly, but in all pairs of sim- 
ulations the Q-parameters are systematically lower after 
TV-body evolution for the runs with feedback compared to 
those without. After lOMyr of TV-body integration this dif- 
ference is most pronounced for Run I (the green lines) , with 
a difference 0.59 in Q-parameter, and least pronounced for 
Run J (the red lines), with a difference 0.02 in Q-parameter 
(though this simulation displays a difference of order 0.7 in 
Q-parameter earlier in the TV-body integration at ~0.5 Myr). 




0.5 



3.3 Stellar evolution 

We show the effects of stellar evolution on our TV-body simu- 
lation^] in Fig. [4] The simulations of clusters with feedback 
are shown by the solid lines, whereas the clusters without 
feedback are shown by the dashed lines. The colours corre- 
spond to the following simulations; red-Run J, green-Run I, 
dark blue-Run UF, cyan-Run UP and magenta-Run UQ. 
In Fig. |4(a)"| we show the evolution of total cluster mass over 



10 



Time (Myr) 



the duration of the TV-body integration, and in Fig. 4(b) 



show the fraction of the initial stellar mass that remains, 
over the 10 Myr of TV-body integration. 

The simulations of clusters that formed with feedback 
have lower initial masses than those that formed without. 
This is because in the simulations without feedback, star for- 
mation is not suppressed by feedback mechanisms, and this 
can lead to the growth of relatively massive stars, which in 
turn r esults in a top-heavy I MF in the clusters without feed- 
back dDale et al ] |2012tl2013t ). However, the formation of the 
more massive stars in the clusters without feedback means 
that a higher fraction of stars lose a significant amount of 
mass (and go supernova) within the duration of the TV-body 
integration. This is most apparent in Runs UF and UQ (the 
dashed dark blue and magenta lines in Fig. |4(b)| respec- 
tively), which both lose more than 50 per cent of their ini- 
tial mass. In contrast, the clusters that formed with feedback 
(the solid lines in Fig. |4(b)[ ) do not contain as many massive 
stars, and typically lose only ~ 20 per cent of their original 



3.4 Half-mass radii 

In Fig. O we compare the evolution of the half-mass 
radii in the simulation pairs. We define the half-mass ra- 
dius based on the centre of mass for the cluster, which 
is problematic if the cluster is in itially substructured (e.g. 
iParker. Goodwin fc Allison! 120111 ). As most of the clusters 
eventually attain a centrally concentrated morphology we 
will examine the differences between simulations, but with 
this caveat in mind. (A better measure of how much the 



2 We also performed the same suite of TV-body simulations with 
stellar evolution switched off, and the results are qualitatively 
similar (such as the evolution of the Q-parameter, evolution of 
the half-mass radius and evolution of surface density). 



Figure 3; Evolu t ion of the Q-parameter 

l|Cartwright fc Whitworthl \2004 ) for all simulations. The 
boundary between centrally concentrated, radially smooth 
clusters (Q > 0.8) and substructured clusters (Q < 0.8) is shown 
by the dotted line. Simulations with feedback are shown by the 
solid lines; simulations with no feedback are shown by the dashed 
lines. The red lines are Runs J, green lines are Runs I, dark blue 
lines are Runs UF, cyan lines are Runs UP and the magenta lines 
are Runs UQ. Runs with feedback preserve structure for longer 
as the cluster evolves, and three simulations with feedback (I, 
UF and UQ) remain substructured beyond 5 Myr. 



cluster expands is perhaps to compare the cumulative dis- 
tributions of stellar s urface densities, as the y do not depend 
on morphology (e.g. IParker fc Meverl [2012). which we will 
explore in Section f3. 61 ) 

There is no obvious pattern in behaviour of the half- 
mass radii. In Run J (the red lines), the feedback-influenced 
cloud loses a substantial (~ 20 per cent) fraction of its gas 
to feedback before the assumed supernova detonation (and 
the subsequent TV-body evolution) , but the resulting cluster 
nevertheless expands more slowly than does the one from 
the corresponding run without feedback. The initial half- 
mass radii for this simulation pair are 1.69 pc for the clus- 
ter that formed with feedback, and 1.20 pc for the cluster 
that formed without. However, after 10 Myr the cluster that 
formed with feedback has a half-mass radius of 5.3 pc com- 
pared to 15.8 pc for the cluster that formed without. A sim- 
ilar result occurs for Run UP (the cyan lines); the initial 
half-mass radii are 3.85 pc (for the cluster that formed with 
feedback) and 3.87 pc (for the cluster that formed without 
feedback). The cluster with feedback (the solid cyan line) 
then expands more slowly and has a final half-mass radius 
of 12.4 pc compared to 19.6 pc for the cluster without feed- 
back (the dashed cyan line). 

For Run I, the cluster that formed with feedback has a 
much larger half-mass radius (3.59 pc) than the simulation 
of the cluster that formed without feedback (0.70 pc) and 
both clusters expand at a similar rate (the green lines in 
Fig-©. However, the high degree of substructure present in 
the simulation with feedback suggests that the determina- 
tion of half-mass radius (which assumes a symmetric, cen- 
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Figure 4. The effects of stellar evolution on the total mass within each simulation. In panel (a) we show the evolution of the total cluster 
mass and in panel (b) we show the fraction of the initial mass of the cluster which is lost through stellar evolution (we include all stars, 
irrespective of whether or not they are still bound). Simulations with feedback are shown by the solid lines; simulations with no feedback 
are shown by the dashed lines. The red lines are Runs J, green lines are Runs I, dark blue lines are Runs UF, cyan lines are Runs UP 
and the magenta lines are Runs UQ. The simulations with no feedback have a more top-heavy IMF, and lose more mass through stellar 
evolution (more than 50 per cent of their initial mass in simulations UF and UQ). 



trally concentrated morphology) could be misleading. Sim- 
ilarly, in Runs UF and UQ (the dark blue and magenta 
lines) the simulations that formed with feedback are highly 
substructured and the small numbers of stars in these simu- 
lations results in the evolution of the half-mass radius being 
rather noisy. This is in part due to the occurrence of su- 
pernovae (compare the dips in half-mass radius with the 



mass-loss in Fig. 4(b) I, which redefines the 'half-mass' on 



very short timescales. In contrast to Run I and Run UP, the 
half-mass radii for the 3 other simulations is lower for those 
clusters that formed without feedback, compared to those 
that formed with feedback. 



3.5 Central densities 

Although ill-defined for substructured clusters, we use the 
half-mass radius to determine a central density, i.e. the stel- 
lar mass density within the half-mass radius. In Fig. HJ] we 
plot the evolution of this density for all five pairs of simu- 
lations. The evolution of simulations that formed with feed- 
back is shown by the solid lines, and the evolution for clus- 
ters that formed without feedback is shown by the dashed 
lines. It is apparent that the initial densit y is lower for clus- 
ters that formed with feedback in the iDale et all i|2012l . 

simulations, than for those that formed without. In 
the most extreme case, Run I that formed with feedback 
(the solid green line in Fig. [5} has an initial density of 
2.1Mqpc -3 , compared to the cluster that formed without 
feedback (the dashed green line), which has an initial den- 
sity of 447 Mq pc -3 . Run UF (the blue lines) also exhibits a 
large difference, with an initial density of 0.3 Mq pc -3 (with 
feedback) compared to 49.9 Mq pc -3 (without feedback). 




Time (Myr) 

Figure 5. Evolution of the half-mass radius for all simulations. 
Simulations with feedback are shown by the solid lines; simula- 
tions with no feedback are shown by the dashed lines. The red 
lines are Runs J, green lines are Runs I, dark blue lines are Runs 
UF, cyan lines are Runs UP and the magenta lines are Runs UQ. 
There is no clear dependence of the evolution of half-mass radius 
on initial conditions (i.e. feedback versus no feedback). 



The subsequent evolution of the densities is non- 
uniform. The clusters that formed without feedback tend 
to become less dense on faster timescales than those that 
formed with feedback, due to rapid mass-loss from stel- 
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Figure 6. Evolution of the central stellar density (i.e. the den- 
sity within the half-mass radius) for all simulations. Simulations 
with feedback are shown by the solid lines; simulations with no 
feedback are shown by the dashed lines. The red lines are Runs 
J, green lines are Runs I, dark blue lines are Runs UF, cyan lines 
are Runs UP and the magenta lines are Runs UQ. Simulations 
that formed with feedback all have lower initial densities than 
those that form without feedback, but the subsequent evolution 
is non-uniform. 



lar winds and supernovae in the iV-body evolution; this is 
most apparent for Runs UF and UQ (the dashed blue and 
magenta lines), which also lose the most stellar mass over 



10 Myr (see Fig. 4(b) I 



3.6 Surface densities 



As discussed in Sections 13.41 and 13.51 determining the half- 
mass radius in a substructured cluster is problematic, and 
a more informative measure of cluster expansion may be 
the evolution of local stellar surface densi ty. We define sur- 
face d ensity, E, based on the definition in lCasertano fc Hutl 
<| 19851 ): 



N ■ 



■kD% 



(2) 



where N is the N th nearest neighbour and Djv is the pro- 
jected distance to that nearest neighbour. We adopt N — 7, 
in order to detect localised pockets of high stellar density but 
not be biased by binary or higher order multiple systems. 

In Fig. [7] we show plots depicting the cumulative dis- 
tributions of surface densities for stars at 0, 1, 5 and 10 
Myr (the solid, dashed, dot-dashed and dotted lines, re- 
spectively) in the iV-body integration of simulation Run J. 
In Fig. 7(a) we show the evolution of the E distribution 



for the cluster that formed with feedback, and in Fig. 7(b) 
we show the evolution for the cluster that formed without 
feedback. The cluster that forms with feedback has an ini- 
tial median surface density E = 141 stars pc -2 , which de- 
creases to 32 stars pc -2 after 1 Myr, 7 stars pc -2 after 5 Myr 
and 2 stars pc -2 after 10 Myr. In contrast, the cluster that 
formed with no feedback has a much higher initial surface 



density, 4518 stars pc~ 2 , which decreases to 26 stars pc -2 af- 
ter 1 Myr, 2 stars pc -2 after 5 Myr and 0.4 stars pc~ 2 after 
10 Myr. The clusters start with very different surface den- 
sities, but converge to similar distributions over the 10 Myr 
of iV-body evolution. 

From inspection of the evolution of the E distribution 
for Run J, it appears that the cluster that formed without 
feedback has been able to form with a much higher surface 
density than the simulation with feedback, as suggested by 
the volume density measurements in Section 13.51 However, 
in Fig. [8] we show the initial surface densities for the simu- 
lation Run I with feedback (panel a) and without feedback 
(panel b) by the solid lines. Recall that the volume density 
for this cluster was 2.1 Mq pc -3 for the cluster that formed 
with feedback, compared to 447 Mq pc -3 for the simulation 
that formed without. However, the median surface densities 
are very similar; the simulation that formed with feedback 
has E = 83 stars pc -2 whereas the simulation that formed 
without has E = 102 stars pc -2 . This is due to the initially 
high level of substructure (Q — 0.38) in the simulation, 
which presents itself as high initial local surface densities 
(the solid line in Fig. 8(a) I, but as a low volume density 



because the volume density is defined from the ill-defined 
cluster 'centre'. 

In Fig. 8(a) we show the evolution of the E distri- 
bution for the cluster that formed with feedback, and in 
Fig. |8(b)| we show the evolution for the cluster that formed 
without feedback. The cluster that forms with feedback 
has an initial median surface density E = 83 stars pc -2 
(the solid line), which decreases to 4 stars pc -2 after 1 Myr, 
0.4starspc -2 after 5Myr and O.lstarspc -2 after 10 Myr 
(the dashed, dot-dashed and dotted lines, respectively). The 
cluster that forms without feedback has a similar initial me- 
dian stellar surface density of 102 stars pc~ 2 , which decreases 
to 10 stars pc -2 after lMyr, 0.8 stars pc -2 after 5 Myr and 
0.3 stars pc -2 after 10 Myr. 

The remaining three pairs of simulations also do not 
show any clear dependence of the evolution of stellar sur- 
face density on whether the clusters formed with or without 
feedback. In all pairs of simulations, the clusters that formed 
without feedback have higher surface densities (although as 
we have seen in Fig. [H] the difference can be marginal), sug- 
gesting that the presence of feedback decreases the density 
of star formation. However, in all simulations, the median 
surface density decreases significantly (often by 3 - 4 orders 
of magnitude) over 10 My r of iV-body in t egrat ion. Such be- 
haviour was also found bv lMoeckel et all (|2012T ) when exam- 
ining the subsequent dynamical evo l ution of SPH simulation 
of star formation bv lBonnell et al.1 (|2008l ). Furthermore, in 
iV-body models of cl usters containing several hundred stars 
in virial equilibrium, [Gieles. Moeckel fc Clarkel l|2012f ) and 
iParker fc Meverl <|2012ri showed that the median surface den- 
sity decreases simply due to two-body relaxation. 



3.7 Bound/unbound stars 

We track the number of stars that remain bound as a func- 
tion of time by calculating the kinetic and potential energies 
for each star. The potential energy of an individual star, Vi, 
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Figure 7. Evolution of stellar surface density for Run J (a) with feedback and (b) no feedback. The cumulative distributions of 
surface density at OMyr (the solid line), 1 Myr (the dashed line), 5Myr (the dot— dashed line) and lOMyr (the dotted line) are shown. 
Before dynamical evolution, the simulation with feedback (a) has a median stellar surface density of 141 stars pc -2 , which decreases 
to 32 stars pc — 2 after 1 Myr, 7starspc — 2 after 5 Myr and 2starspc~ 2 after 10 Myr. In contrast, the simulation with feedback (b) has 
an initial median stellar surface density of 4518 stars pc -2 , which decreases to 26 stars pc -2 after 1 Myr, 2 stars pc -2 after 5 Myr and 
0.4 stars pc~ 2 after 10 Myr. 




Figure 8. Evolution of stellar surface density for Run I (a) with feedback and (b) no feedback. The cumulative distributions of surface 
density at OMyr (the solid line), 1 Myr (the dashed line), 5 Myr (the dot-dashed line) and 10 Myr (the dotted line) are shown. Before 
dynamical evolution, the simulation with feedback (a) has a median stellar surface density of 83 stars pc — 2 , which decreases to 4 stars pc — 2 
after lMyr, 0.4starspc — 2 after 5Myr and O.lstarspc - 2 after lOMyr. The simulation with feedback (b) has a similar initial median 
stellar surface density of 102 stars pc -2 , which decreases to 10 stars pc — 2 after 1 Myr, 0.8 stars pc — 2 after 5 Myr and 0.3 stars pc -2 after 
10 Myr. 



is given by: 

(3) 

■ / ■ ^ij 

where rtii and rrij are the masses of two stars and is the 
distance between them. The kinetic energy of a star, Ti is 



given thus: 

1 2 
Ti = -TOi|vi ~v c i| , (4) 

where and v c i are the velocity vectors of the star and the 
centre of mass of the cluster, respectively. A star is bound 
if Ti + Vi < 0. 
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Figure 9. Evolution of the number fraction of bound stars. Sim- 
ulations with feedback are shown by the solid lines; simulations 
with no feedback are shown by the dashed lines. The red lines are 
Runs J, green lines are Runs I, dark blue lines are Runs UF, cyan 
lines are Runs UP and the magenta lines are Runs UQ. There is 
no correlation between the evolution fraction of bound stars and 
the initial conditions (i.e. feedback versus no feedback). 



In Fig. [9] we show the number fraction of stars that re- 
main bound over the iV-body integration for all five pairs of 
simulations. Again, the simulations that formed with feed- 
back are shown by the solid lines, and the simulations that 
formed without feedback are shown by the dashed lines. 
The colours correspond to the following simulations; red- 
Run J, green-Run I, dark blue-Run UF, cyan-Run UP and 
magenta-Run UQ. 

Firstly, we note that the fraction of stars that are bound 
at the start of the iV-body integration is less than 1 in all 
of the simulations, even those that are initially virialised, or 
slightly subvirial. This is due to the initially non-spherical, 
substructured nature of the simulations, which makes the 
definition of the centre of mass of the cluster problematic, 
and also that the individual velocities of stars which form 
on the outskirts of the cluster can be high. 

The simulations that formed from unbound clouds (UF, 
UP and UQ) have initially lower fractions of bound stars 
than the simulations that formed from bound clouds (Runs J 
and I). 

The fraction of stars that are still bound at the instant 
of the removal of remaining gas appears to be slightly depen- 
dent on whether the cluster forms with or without feedback. 
In all of the simulation pairs, the cluster that formed with 
feedback has a lower fraction of bound stars initially than the 
cluster that formed without feedback (the solid and dashed 
lines in Fig. [5] respectively). However, in Runs J, I and UQ 
(the red, green and cyan lines) the cluster that formed with 
feedback retains more bound stars over the 10 Myr of N- 
body integration. The remaining two simulations (UF and 
UQ) have similar final fractions of bound stars, irrespective 
of whether the clusters formed with or without feedback. 



3.8 Mass segregation 

We use two different algorithms to search for mass segrega- 
tion in the pairs of simulations. Firstl y, we apply the Min i- 
mum Spanning Tree (MST) method of lAllison etall l|2009l) . 
which compares the MSTs of random sets of stars with 
the MST of a chosen subset. The ratio of the average ran- 
dom MST to the subset, Amsr, is significantly larger than 
unity if mass segregation is present (such as is observed in 
the Orion Nebula Cluster, iHillenbrand fc Hartmannlll998l ; 
lAllison et al.ll2009l ). or less than unity if the most massive 
stars are more spread out than the average stars. The sec- 
ond al gorithm is the m — E method of lMaschberger fe Clarke] 
l|2011f ). which plots the surface density, E, of stars against 
their mass, m. Mass segregation is defined as the most mas- 
sive stars having significantly larger local stellar surface den- 
sities compared to the average stars in the cluster. 

There are two approaches to measuring mass segrega- 
tion in clusters. One can either measure the global level 
of segregation (i.e. including all stars in the cluster - e.g . 
lAllison et ai1l2009l ; lOlczak et ai1l201ll ; iParker et all l201lK . 
or alternatively, identify subclusters and measure mass 
segregation on localised scales (e.g. Kirk fc Mversl l201ll ; 
iMaschberger fc Clarke! 1201 ll ; iGirichidis et al.ll2012l ). Here, 
we focus on global mass segregation, because we expect the 
subsequent iV-body evolution to facilitate two-body relax- 
ation and potentially lead to dynamical mass segregation 
when the stars are well-mixed. 

We search for mass segregation at Myr, 1 Myr and 
5 Myr in the simulations. Several clusters display hints of 
mass segregation, but not on a significant level. We there- 
fore conclude that none of the clusters exhibit primordial, 
or subsequent dynamical mass segregation. 



4 DISCUSSION 

The main result in Section [3] can be summarised as follows. 
The evolution of cluster structure, as determined by the Q- 
parameter, demonstrates that clusters which form under the 
influence of feedback retain primordial structure for longer 
th an clusters that form w ithout feedback in the simulations 
of iDale et ail (|2012l . l2013l ). Furthermore, three clusters that 
form with feedback remain substructured throughout the 
subsequent iV-body evolution. This result in itself is not 
particularly helpful, as the differences in the amount of sub- 
structure can be merely relative; in two pairs of simulations 
(Runs J and UP) both clusters attain a centrally concen- 
trated morphology, but the clusters that formed with feed- 
back have lower Q-parameters. 

Furthermore, the evolution of the Q -parameter can be 
highly stochastic. I Parker fc Meverl )|2012l ) demonstrate that 
an initially highly substructured (Q < 0.5), subvirial cluster 
(T/|l/| = 0.3, where T and V are the total kinetic and po- 
tential energies of the stars in a cluster) erases its substruc- 
ture within 0.5 Myr as it becomes centrall y concentrated due 
to violent relaxation (|Allison et al.ll2010l ). An initially sub- 
structured and supervirial cluster (T/|V| = 1.5) expands 
and generally retains its substructure, although 1 in 10 clus- 
ters (with statistically identical initial conditions) become 
centrally concentrate d due to the stochast ic nature of the 
dynamical evolution (|Parker fc Meverl [2012). In both evolu- 
tionary scenarios, the range in final Q-parameters between 
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10 identical simulations is ~ 0.5, implying that the differ- 
ences between the simulations with and without feedback 
presented here could just be due to the stochastic dynami- 
cal evolution of the clusters over 10 Myr. 

Howe ver, there are several diffe rences between the sim - 
ulations in I Allison et all l|2010l ) and lParker fc Meveij (|2012l ) 
and the clusters we dynamically evolve here. Firstly, the 
fractal simulations in those pape rs assume highly cor- 
relate d velocities on local scales (|Goodwin fc Whitworthl 
l2004h . which increases the amount of dynamical mixing that 
takes place as the substructure violently relaxes. Secondly, 
those simulations were of fractal clusters with N = 1000 
stars and radii of ~lpc, (with a typical central density of 
~ 2000 Mq pc -3 following the subviri al collapse) wherea s the 
most compact SP H simulation from iDale et al l <|2012l ) and 
IDale et all (|2013h is the cluster in Run J that formed with- 
out feedback, which formed 578 stars and has a radius of 
~lpc (and central density of 220 Mq pc -3 ). At the other 
end of the scale, Run UQ formed only 80 stars and has a 
radius of ~12pc (and a central density of 0.6 Mq pc -3 ). 

Therefore, we expect the crossing time (and the related 
relaxation time) to be much longer for the simulations pre- 
sented here. As an example, we co mpute the relaxation time, 
^rciax, of a typical simulati on from I Allison et al.l (I2010T ). us- 
ing the following formula (|Binnev fc Tremainell 19871 ): 

N 

^rclax = 7T"j TT^cross, (5) 

o In iV 

where N is the number of stars an d Across is th e cross - 
ing time. For a typical simulation in lAllison et all |2010l ). 
Across = 0.1 Myr and N = 1000, which corresponds to a 
relaxation time t rc iax = 2 Myr. On the other hand, the 
crossing time for the cluster in Run J that formed with 
feedback is t CIOSS ~ 0.7 Myr, and its relaxation time (for 
N = 685 stars) is 9 Myr. Conversely, the cluster in Run J 
that formed without feedback has a smaller radius (but sim- 
ilar velocity dispersion) and hence a much shorter crossing 
time (icross ~ 0.25 Myr), and a relaxation time of 2.8 Myr. 
This is evident in the behaviour of the Q-parameter; in 
Fig. [3] the cluster in Run J that formed without feedback 
(the red dashed line) has a constant Q-parameter after 
~1.5 Myr, whereas the simulation that formed with feedback 
(the red solid line) has a Q-parameter that is still increasing 
at 10 Myr. 

A similar dependence of the evolution of the Q- 
parameter on the relaxation time is found for each pair of 
simulations; in all cases, the clusters that formed with feed- 
back have longer relaxation times, which is likely due to the 
initial densities of the clusters. In Fig. [6] we show that the 
initial volume density is lower for the simulations that form 
with feedback, due to gas expulsion, regulation of stellar 
mass growth, and the triggering and redistribution of star 
formation. We suggest that the two clusters that form with 
feedback which never erase their substructure are simply too 
diffuse to do so (Runs UF and UQ, which have initial densi- 
ties less than 1 Mq pc -3 ). This is also echoed by the initial 
surface densities of these clusters, which are 0.6 stars pc -2 
for Run UF with feedback and 2 stars pc -2 for Run UQ with 
feedback. It is not possible to infer the initial density after 
1 Myr of dynamical evolution as the clusters all expand and 
relax, and clusters with very different initial densities evolve 
to have very similar densities (Figs. [6] [7]and[8]|. 



Whilst it is possible to infer differences between the sim- 
ulations with and without feedback, one needs to make an 
estimate of the initial density in order to calibrate the evo- 
lution of the Q-parameter. As we have seen, a Q-parameter 
which indicates a centrally concentrated morphology at an 
age of > 2 Myr could have formed with or without feedback. 
The only constraint we can place is that all the cl usters that 
formed without feedback in the simulations of IDale et al.l 
(2012) and Dal e et ail (|2013l ) become centrally concentrated. 
If a cluster is observed at ^5 Myr to be substruc tured, it 
is either highly supervirial (jParker fc Meverll2012i ). and/or 
formed under the influence of feedback. 

Other diagnostics of dynamical evolution do not distin- 
guish between feedback and non-feedback star formation. 
Initial differences in the IMF (which is more top-heavy for 
simulations which are not regulated by feedback) are re- 
duced by subsequent stellar evolution. Additionally, it is 
not possible to distinguish the two types of star formation 
by considering the kinematics of the clusters' evolution, as 
there is no systematic dependence on the fraction of bound 
stars. 



4.1 Scaled velocities 

In addition to our default initial conditions, we also per- 
formed a set of simulations in which we scaled the initial 
bulk motion of the stars in each cluster to be in virial equi- 
librium. Qualitatively, the results do not depend on the ini- 
tial velocity set-up. There is again no trend in the evolu- 
tion of the half-mass radius, bound fraction, surface density, 
depending on whether the cluster formed with or without 
feedback. We show the evolution of structure, as measured 
by the Q-parameter, for these (initially) virialised clusters 
in Fig. 1101 As in the case of the simulations with non-scaled 
velocities, the clusters that form with feedback retain more 
structure during the subsequent dynamical evolution. 

Comparison of Fig.[TU]to Fig.[3]shows that the different 
velocity scaling has not made much difference to the evolu- 
tion of the Q-parameter. The simulation pair with the most 
stars (Run J - the red lines) does not become as centrally 
concentrated when in virial equilibrium as the unsealed sim- 
ulations in Fig. [3] which may be due to the slightly subvirial 
initial velocities in the non-scaled s imulations (as subviria l 
clusters reach higher 0-parameters. [Parker fc Meverir2012l ). 
The remaining clusters follow similar evolution to those in 
the simulations that have the stellar velocities scaled to be in 
virial equilibrium, apart from Run I with feedback (the solid 
green line) which reaches the border between structured and 
centrally concentrated earlier when in virial equilibrium (at 
5 Myr instead of 8 Myr) . 

All other parameters (density evolution, fraction of 
bound stars) are similarly unchanged when the velocities 
of stars are scaled to be in virial equilibrium. 



5 CONCLUSIONS 

We have tak en the end states of five pairs of hydrodynamical 
simulations l|Dale et al.ll2012j.l2013h and evolved them as N- 
body simulations for 10 Myr to look for differences between 
the clusters in the simulations which formed under the in- 
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Figure 10. Evolution of the Q-parameter for all simulations in 
which we scaled the velocities to be initially in virial equilibrium. 
The boundary between centrally concentrated, radially smooth 
clusters (Q > 0.8) and substructured clusters (Q < 0.8) is shown 
by the dotted line. Simulations with feedback are shown by the 
solid lines; simulations with no feedback are shown by the dashed 
lines. The red lines are Runs J, green lines are Runs I, dark blue 
lines are Runs UF, cyan lines are Runs UP and the magenta 
lines are Runs UQ. The trend in the evolution of substructure 
is similar to the simulations with non-scaled velocities; however, 
the spread in Q between simulations is lower, and Run I with 
feedback (the solid green line) loses its substructure earlier than 
in the non-scaled velocity case (at 5 Myr instead of 8 Myr). 



whilst none of the non-feedback simulations are able to sup- 
port substructure for more than 1 Myr of iV-body evolution. 

This difference is a straightforward consequence of the 
feedback clusters having lower volume densities and conse- 
quently longer relaxation times, so that in simulations UF 
and UQ, there is in fact little stellar mixing at all. These 
lower volume densities are in turn due to prolonged grad- 
ual mass loss from the parent clouds, and to the triggering 
and redistribution of star formation outside the principal 
subclusters. 

However, the feedback and non-feedback clusters have 
varying degrees of initial substructure, and there is no clear 
trend for the feedback calculations to have more or less 
initial substructure than their counterparts. It is thus very 
difficult to draw reliable conclusions about the influence 
feedback has had on a particular system from the observed 
Q-parameter at any epoch. The strongest statement that 
could be made is that a substructured cluster > 5 Myr in 
age is likely to have experienced strong feedback. However, 
even this inference requires knowledge of the initial density 
of the cluster, and its initial virial state, and is based on a 
small number of simulations of small- N clusters. 
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fluence and feedback, compared to the evolution of clusters 
that formed without feedback. 

Gas expulsion and the other effects of photionization 
feedback - the triggering, aborting and redistribution of 
star formation - are expected to affect the initial dynamical 
states of young clusters, since they influence the masses, po- 
sitions and velocities of the stars. The five pairs of iV-body 
simulations presented here demonstrate clearly that this is 
the case, by comparing the subsequent evolution of clusters 
which either have or have not experienced photoionization 
feedback before supernovae are supposed to instantaneously 
clear them of gas. 

Although the evolution of the feedback and no-feedback 
members of each pair of simulations is different, it is dif- 
ficult to identify general trends in these differences - for 
example, some feedback-influenced clusters lose stars more 
rapidly than their counterparts during their iV-body evolu- 
tion and some do not. 

The only general result that we have been able to 
identify is that the feedback-i nfluenced cl usters which 
formed in the simulations of bale et alj (|2012l . I2013T ) 
tend to lose their substructure (as measu r ed by the 
Q-pa rameter, ICartwright fc Whitworthl |2004| ; ICartwrightl 
2009) more slowly than their non-feedback-influenced op- 
posite numbers. The most strongly feedback-aff e cted simu- 
lations - Runs I , UQ and UF from iDale et al.l (|2012l ) and 
iDale et al l |2013l ) - retain some substructure for 5-10 Myr, 
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